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We have investigated electronic and magnetic properties of graphene nanodisks (nanosize triangular 
graphene) as well as electromechanical properties of graphene nanoj unctions. Nanodisks are nanomagnets 
made of graphene, which are robust against perturbation such as impurities and lattice defects, where the ferro- 
magnetic order is assured by Lieb's theorem. We can generate a spin current by spin filter, and manipulate it by 
a spin valve, a spin switch and other spintronic devices made of graphene nanodisks. We have analyzed nan- 
odisk arrays, which have multi-degenerate perfect flat bands and are ferromagnet. By connecting two triangular 
graphene corners, we propose a nanomechanical switch and a rotator, which can detect a tiny angle rotation by 
measuring currents between the two corners. By making use of the strain induced Peierls transition of zigzag 
nanoribbons, we also propose a nanomechanical stretch sensor, in which the conductance can be switch off by 
a nanometer scale stretching. 
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I. INTRODUCTION 

Graphene, which is a one-layer thick honeycomb structure 
of carbon, is an amazing material^ . Electrons exhibit high mo- 
bility and travel micron distances without scattering at room 
temperature. It is a very thin and strong material, showing 
very high thermal conductivity. Graphene is now a main topic 
of nanoscience. 

Much attention has been focused on graphene nanorib- 
bons, which is a one-dimensional ribbon-like derivatives of 
graphene. They have various band structure depending on the 
edge and width. In particular, zigzag gaphene nanoribbons 
show edge ferromagnetism due to almost flat low-energy band 
at the Fermi level. There are a profusion of papers on them, 
among which we cite some of early works^Tii. Another ba- 
sic element of graphene derivatives is a graphene nanodisk-. 
It is a nanometer- scale disk-like material which has a closed 
edge. It may be considered as a giant molecule made of 
aromatic compound. It is possible to manufacture them by 
etching a graphene sheet by Ni nanoparticles^. Among them, 
trigonal zigzag nanodisks have a novel electric property that 
there exist half-filled zero-energy states in the non-interacting 
regime, as was revealed first by the tight-binding model^ 
and then by first-principle calculations^"^. Various remark- 
able properties of nanodisks have been investigated exten- 
sively in a series of works^^i^—. Nanodisk is also referred to 
as nanoisland^, nanoflake^iJ^ii^, nanofragmenU^ or graphene 
quantum dot^^'^^. 

Nanoribbons and nanodisks correspond to quantum wires 
and quantum dots, respectively. They are candidates of fu- 
ture carbon-based nanoelectronics and spintronics alternative 
to silicon devices. A nanoribbon-nanodisk complex can in 
principle be fabricated, embodying various functions, only 
by etching a graphene sheet. Furthermore, graphene is com- 
mon material and ecological. In this paper, exploring elec- 
tronic, magnetic and electromechanical properties of trigonal 
zigzag graphene nanodisks [see Fig.l], we propose some ap- 
plication of nanodisk-nanoribbon complex to nanoelectronics, 
spintronics and electromechanics devices 



(b) nanodisk-lead system 




FIG. 1 . (a) Geometric configuration of a trigonal zigzag nanodisk. 
The nanodisk size is defined hy N — A/ben — 1 with N^en the number 
of benzenes on one side of the trigon. Here, N — 5. The A (B) sites 
on the lattice are indicated by red dots (while circles). The electron 
density is found to be localized along the edges, (b) The nanodisk- 
lead system. A nanodisk is connected to the right and left leads by 
tunneling coupling. It may act as a spin filter. 



II. ENERGY SPECTRUM 

We calculate the energy spectrum of the nanodisk based on 
the nearest-neighbor tight-binding model, which has proved to 
describe accurately the electronic structure of graphene, car- 
bon nanotubes, graphene nanoribbons and other sp^ carbon 
materials. The Hamiltonian is defined by 
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where Si is the site energy, tij is the transfer energy, and c] 
is the creation operator of the tt electron at the site i. The 
summation is taken over all nearest neighboring sites (i^j). 
Owing to their homogeneous geometrical configuration, we 
may take constant values for these energies, Si = £f and Uj = 
t ^ 2.70eV. There exists one electron per one carbon, and the 
band-filling factor is 1/2. Then, the diagonal term yields just 
a constant, Ey^Nq, and can be neglected in the Hamiltonian, 
where A^c is the number of carbon atoms. 

We define the size A" of a nanodisk by A' = Aben — 1, where 
Aben is the number of benzenes on one side of the trigon as in 
Fig. 1(a). It can be shown^ that the determinant associated with 
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the Hamiltonian ([T]) has such a factor as 

det[€l-H{Nc)] (xe^, (2) 

implying A/^-fold degeneracy of the zero-energy states. The 
gap energy is as large as a few eV for nanodisks with small N, 
where it is a good approximation to investigate the electron- 
electron interaction physics only in the zero-energy sector, by 
projecting the system to the subspace made of those zero- 
energy states. As we shall see in SectionjVl the approximation 
remains to be good even for those with large N. 

III. TRIGONAL SYMMETRY 




FIG. 2. Probability density flow for the zero-energy states in the 
nanodisk with N = 7. The representation of the trigonal symmetry 
group C3v is indicated in the parenthesis. A vortex appears at the 
center of mass for the state belongs to the E (doublet) representation. 
The winding number at the center is 2 in the state \ki). 



The symmetry group of a trigonal nanodisk is Csv , which is 
generated by the 27r/3 rotation Cs and the mirror reflection . 
It has the representation {Ai, A2, E}. The Ai representation 
is invariant under the rotation C3 and the mirror reflection . 
The A2 representation is invariant under C3 and antisymmetric 
under at,. The E representation acquires ±27r/3 phase shift 
under the 27r/3 rotation. The Ai and A2 are 1 -dimensional 
representations (singlets) and the £^ is a 2-dimensitional rep- 
resentation (doublet). These properties are summarized in the 
following character table: 
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The zero-energy sector consists of N orthonormal states. 
We are able to index them by the wave number k along 
the edge. It is a continuous parameter for an infinitely long 
graphene edge. According to the tight-binding-model result, 
the flat band emerges for 

- TT < afc < -27r/3 and 27r/3 < aA: < tt. (4) 

We focus on the wave function ijj {x^y) at one of the A sites 
on an edge, and investigate the phase shift when we step over 
to the neighboring site [see Fig. 1(a)]. There are N links along 
one edge of the size-TV nanodisk, for which we obtain the 
phase shift Nak, where a is the spacing between the neighbor- 
ing A sites and k is the wave number along the edge. On the 
other hand, the phase shift is tt at the corner. The total phase 
shift is 3Nak + Sir, when we encircle the nanodisk once. By 
requiring the single- valueness of the wave function, it is found 
to be quantized as 

akn = ± [(2n + 1)/37V + 2/3] tt, (5) 

where it follows that < n < {N — l)/2 from the allowed 
region of the wave number. 

We may group the states according to the trigonal symme- 
try ([3]). With respect to the rotation there are three elements 
C3, C3, C3, which correspond to 1, e^^*/^, e^^*/^. Accordingly, 
the phase shift of one edge is 0, 27r/3, 47r/3. The state is 



grouped according to the representation of the trigonal sym- 
metry group Csy as follows. 
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A2 (singlet): - | - fc^), ^ 

E (doublet): \-k^), for k^ 
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The zero-energy state is indexed by the quantized wave num- 
ber as \k^) with ©. 

To see the meaning of the wave number k!^ more in detail^, 
we have calculated the probability density flow. 
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for states which we show for the case of = 7 in Fig. 2. 
We observe clearly a texture of vortices. These vortices mani- 
fest themselves as magnetic vortices perpendicular to the nan- 
odisk plane when the electromagnetic fields are coupled. 
The total winding number Ayortex is calculated by 
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with m = 0, 1, 2, • • • , [{N - l)/2j in the size-A^ nanodisk, 
where [aj denotes the maximum integer equal to or smaller 
than a. We find Ayortex = 3n for /cJJ, Ayortex = 3n + 1 for 
k~j^-^ and Ayortex = 3n + 2 for The wave functions are 
classified in terms of modulo of the total winding number: The 
wave function belongs to the E-representation and has chiral 
edge mode for Avortex = 1,2 (mod 3), and belongs to the A- 
representation and has non-chiral edge mode for AVortex = 
(mod 3). The winding number of the vortex at the center of the 
nanodisk is 0, 1, 2 in the state respectively. 



IV. QUASIFERROMAGNET 

The total spin of the ground state is determined by Lieb's 
theorem. The total spin is given by the difference of the 
A site and the B site, S = ^\Na- Nb\, where Na and 
Nb are number of A site and B site. Here, S = ^A" in 
the size- A' nanodisk, where Na = {N ^ 1){N + 6)/2 and 
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FIG. 3. Thermodynamical properties of the nanodisk-spin system, 
(a) The specific heat C in unit of k^N. (b) The susceptibihty x in 



unit of Sg. The size is = 1, 2, 2^ 



The horizontal axis 



stands for the temperature T in unit of JN/ . The arrow represents 
the phase transition point Tc in the Umit N ^ oo. 



Nb = {N^ 2) {N + 3)/2. Hence we expect a nanodisk to act 
as a ferromagnet. The ferromagnetic ground state is robust 
against perturbations such as randomness and lattice defects 
since it is assured by Lieb's theorem. This feature brings out a 
remarkable contrast between nanodisks and nanoribbons. The 
ferromagnetic order is fragile due to the lack of Lieb's theo- 
rem in the case of nanoribbons. 

We investigate ferromagnetic properties by introducing 
Coulomb interactions into the zero-energy sector^. We have 
calculated specific heats and susceptibilites at temperature 
T in Fig. 3. There appear singularities in thermodynamical 
quantities sls N oo, which represent a phase transition 
at Tc between the ferromagnet and paramagnet states, where 
Tc = JN/2kB. For finite N, there are steep changes around 
Tc, though they are not singularities. It is not a phase tran- 
sition. However, it would be reasonable to call it a quasi- 
phase transition between the quasiferromagnet and paramag- 
net states. Such a quasi-phase transition is manifest even in 
finite systems with TV = 100 1000. 

The specific heat takes nonzero- value for T > Tc, as shown 
in Fig. 3 (a), which is zero in the limit N ^ oo. The result in- 
dicates the existence of some correlations in the paramagnet 
state. On the other hand, the susceptibility x always shows 
the Curie- Weiss low x 1/^ near T = 0, and exhibits also 
a behavior showing a quasi-phase transition at T = Tc, as 
shown in Fig. 3(b). In the finite system, the expectation value 
of Sz,tot is always zero because there is no spontaneous sym- 
metry breakdown in the finite system, and the behavior is that 
of paramagnet. 



V. MAGNETISM OF LARGE NANODISKS 

For large N nanodisks, the band gap decreases inversely 
proportional to the size. One may wonder if our analysis based 
only on the zero-energy sector is valid. Indeed, the size of 
experimentally available nanodisks is as large as A/" = 100 ~ 
1000. We wish to argue that our analysis based on the zero- 
energy sector is essentially correct, even if the size N of the 
nanodisk is large and the band gap becomes very narrow. 

Near the Fermi energy, the density of states (DOS) D (e) 
consists of that of the bulk graphene and an additional peak 
at the zero-energy states due to the edge states for 1, 




FIG. 4. The density of states of (a) graphene, (b) graphene nanorib- 
bon and (c) graphene nanodisk. The horizontal axis is the energy e 
in unit of t. 



as illustrated in Fig. 4. Hence, together with spin degrees of 
freedom, it behaves as 



D{e) = 2cNc\e\^2NS (e) , 



(9) 



with a certain constant factor c. The linear term is due to the 
bulk states, and the Dirac delta function term is due to the 
edge states. The important point is that the edge- state peak 
is clearly distinguished from the DOS due to the bulk part. It 
is enough to take into account only the zero-energy sector to 
analyze physics near the Fermi energy, since the contribution 
from the edge states is dominant. 

We calculate the magnetization of a nanodisk when its size 
is large. We start with the Hubbard Hamiltonian, 

H = ^e{k) cl^Cka + ^ ^ ^UqA^-qi^k'iCk^. (10) 
ka kk'q 

Let (n^) , (n^) be the average numbers of the up and down 
spins. The magnetization is given by (m) = (n^) — (n^). It 
is determined self-consistently by the relation 

{m) = JdsD{s)[f{s-A)-f{s + A)], (11) 

in terms of the Fermi distribution function 

/(x) = (exp[(x-M)ABT] + l)-i (12) 

and 

A = JL^m) + lh. (13) 

Substituting the formula ^ into the Storner equation (fTTl) . we 
obtain 



4Li2 (- 



(14) 

with the dilogarithm function Li2 (x). It is difficult to solve 
this equation for (m) self-consistently at general temperature 
T. We examine two limits, T ^ and T ^ oo. 

For the zero temperature (T 0) we obtain the magneti- 
zation as 

(m) = TV + {mf + ^ (m) + 0{h'). (15) 

Because |(m)| < N, it follows that (m) = N ^ 0(1). The 
contribution from the bulk gives a negligible correction to the 



FIG. 5. (a) An electron tunnels from the left lead to the nanodisk 
and then to the right lead. The system is a reminiscence of a metal- 
ferromagnet-metal junction, (b) Only electrons with the same spin 
direction as the nanodisk spin can pass through the nanodisk freely. 
As a result, when we apply a spin-unpolarized current to the nan- 
odisk, the outgoing current is spin polarized to the direction of the 
nanodisk spin. Consequently, this system acts as a spin filter. 



FIG. 6. The spin valve is made of two nanodisks with the same 
size, which are connected to leads. Applying an external magnetic 
field, we control the spin direction of the first nanodisk to be \0) 
and that of the second nanodisk to be |t). The incoming current is 
unpolarized, but the outgoing current is polarized and its magnitude 
can be controlled continuously. This acts as a spin valve. 



total magnetization. Hence the magnetization is (m) = 
and the ground state is fully poralized whenever U 7^ 0. 
Ferromagnetism occurs irrespective of the strength of the 
Coulomb interaction. The magnetization is propotional not 
to A^c but N. In this sence the ground state of nanodisk is not 
bulk ferromagnet but surface ferromagnet, which is consistent 
with the previous result. 

We next investigate the high temperatuer limit (T 00). 
Using the Taylor expansion of the dilogarithm function, 



Li2(- 



we find 
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The leading term is the second term, and hence the main con- 
tribution comes from the bulk. The solution is only (m) = 
for which A = 0. There is no magnetization at high tempera- 
ture. 

A comment is in order. We have assumed that the mag- 
netization axis is fixed and only longitudinal fluctuations of 
the magnetic moments take place. In general, spin-wave-like 
fluctuations are dominant at the edges of graphene^^, because 
they are gapless Goldstone modes. On the contrary, there exist 
no gapless Goldstone modes in graphene nanodisks, because 
the edge is finite and closed. Furthermore, when the length 
of the edge is very small, spin-wave-like fluctuations have a 
large gap. Hence, our approximation is valid for nanodisks. 



VI. APPLICATION FOR SPINTRONIC DEVICES 

The nanodisk-spin system is a quasiferromagnet, which is 
an interpolating system between a single spin and a ferro- 
magnet. It is easy to control a single spin by a tiny current 
but it does not hold the spin direction for a long time. On 
the other hand, a ferromagnet is very stable, but it is hard 
to control the spin direction by a tiny current. A nanodisk 



quasiferromagnet has an intermediate nature: It can be con- 
trolled by a relatively tiny current and yet holds the spin di- 
rection for quite a long time^. Taking advantage of these 
properties we have already proposed elsewhere some appli- 
cations of graphene nanodisk-lead systems to spintronic de- 
vices. They are spin filter, spin memory, spin amplifier, spin 
valve, spin-field-effect transistor, spin diode and spin switch, 
among which here we make a consice review of spin filter, 
spin valve and spin switch. We newly propose nanodisk ar- 
rays and nanomechanical switch. 

Spin filter: We consider a lead-nanodisk-lead system [see 
Fig. 1(b)], where an electron makes a tunnelling from the left 
lead to the nanodisk and then to the right lead. This sys- 
tem is a reminiscence of a metal-ferromagnet-metal junction 
[see Fig. 5]. If electrons in the lead has the same spin di- 
rection as the nanodisk spin, they can pass through the nan- 
odisk freely. However, those with the opposite direction 
feel a large Coulomb barrier and are blocked (Pauli blockade 
[Fig.5(c)])^^ As a result, when we apply a spin-unpolarized 
current to the nanodisk, the outgoing current is spin polarized 
to the direction of the nanodisk spin. Consequently, this sys- 
tem acts as a spin filter. 

Spin valve: A nanodisk can be used as a spin valve, in- 
ducing the giant magnetoresistance effect. We set up a sys- 
tem composed of two nanodisks sequentially connected with 
leads [see Fig. 6]. We apply external magnetic field, and 
control the spin direction of the first nanodisk to be |^) = 
cos I It) + sin I II), and that of the second nanodisk to be 
|0) = It). We inject an unpolarized- spin current to the first 
nanodisk. The spin of the lead between the two nanodisks is 
polarized into the direction of\9). Subsequently the current is 
filtered to the up- spin one by the second nanodisk. The out- 
going current from the second nanodisk is I^^^ = / cos | . We 
can control the magnitude of the up-polarized current from 
to / by rotating the external magnetic field. The system act as 
a spin valve. 

Spin switch: We consider a chain of nanodisks and leads 
connected sequentially [see Fig. 7]. Without external magnetic 
field, nanodisk spins are oriented randomly due to thermal 
fluctuations, and a current cannot go through the chain. How- 
ever, when and only when a uniform magnetic field is applied 
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FIG. 7. A chain of nonodisks and leads acts as a spin switch. With- 
out external magnetic field, nanodisk spins are oriented randomly 
due to thermal fluctuations, and a current cannot go through the 
chain. However, as soon as a uniform magnetic field is applied to 
all nanodisks, the direction of all nanodisk spins become identical 
and a current can go through. 




FIG. 8. (a) Illustration of one-dimensional nanodisk arrays, (b) 
Band structure of one-dimensional nanodisk arrays. (c) Illustration of 
two-dimensional nanodisk arrays. 



to all nanodisks, the direction of all nanodisk spins become 
identical and a current can go through. Thus the system acts 
as a spin switch, showing a giant magnetoresistance effect. 
The advantage of this system is that a detailed control of mag- 
netic field is not necessary in each nanodisk. 

Nanodisk arrays: We investigate nanodisk arrays, which 
are materials where nanodisks are connected in one- or two- 
dimentions. These structures have already been manufactured 
by etching a graphene sheet by Ni nanoparticles^. We show 
an example of a trigonal-zigzag-nanodisk array sharing one 
zigzag edge as in Fig. 8 (a). We show the corresponding band 
structure in Fig. 8(b). It is intriguing that there are A^-fold de- 
generate perfect flat bands in the nanodisk with size N. This 
fact is confirmed by Leib's theorem. Each nanodisk has spin 
N/2 and makes ferromagnetic coupling between two nan- 
odisks. In the same way we can make two-dimensional nan- 
odisk arrays. It is to be emphasized that they show ferromag- 
netism, and not quasiferromagnetism, though they are made 
of nonmagnetic materials. The perfect flat band will be robust 
even when electron interactions are introduced. 

Nanomechanical switch: We construct a nanomechan- 
ical switch contacting two graphene trigonal corners [see 
Fig. 9(a)]. We assume the angle between two corners is 6. This 
angle is tuned by an external mechanical force. The carbon- 
skeleton structure is made of a-bonds, and is very rigid except 
for this rotational degree of freedom. 




On state Off state 




FIG. 9. (a) Illustration of a nanomechanical junction made of two 
triangular graphene corners. There are on and off states, switched by 
a rotation between two cornes. (b) Illustration of a nanomechanical 
rotator made of two corners and a central rhombus rotating rather 
freely. 

The conductance is determined by the overlap integral of 
TT-electrons between two corners, which is given by 

\{Pz COS ^py sinO \pz)\ = cos^ 0. (18) 

When the two planes are parallel (0 = 0), the overlap takes 
the maximum value and 7r-electrons can go through the con- 
tact. This is the on state. When the two planes are orthogo- 
nal (0 = 7r/2), the overlap takes the minimum value and tt- 
electrons can not go through the contact. This is the off state. 
The angle is changed nanomechanically. The system acts as a 
nanomechanical switch. It could detect the angle very sensi- 
tively and be useful for detect nanomechanical oscillations. 

By connecting two nanomechanical junctions, we can con- 
struct a nanomechanical rotator [Fig. 9(b)], where two corners 
are suspended mechanically while the central rhombus rotates 
rather freely. This structure may be useful to detect molecular 
dynamics. When molecules contact the rotator, they are de- 
tected by rotating it and changing the resistance between the 
two corners. 

Peierls instability: We connect two triangular corners 
with a zigzag nanoribbon [see Fig. 10(a)]. This structure 
has already been manufactured experimentally^. Poly acety- 
lene has the Peierls instability: Carbons with conjugate 
bonds are spontaneusly deformed into alternating single and 
double bonds. We expect the Peierls instability to occur 
also in graphene nanoribbons [Fig. 10(b)], because grapahene 
nanoribbon is a natural extension of poly acetylene^. (We de- 
note the width of a nanoribbon by W with 1^ = for poly- 
acetylene and = 1 for polyacene.) However, as we now 
show, the Peierls instability will not occur in nanoribbons with 
the width 1^ > 2. On the contrary, it is possible to induce the 
Peierls transition manually by streaching a nanoribbon. Mak- 
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FIG. 10. (a) Illustration of a configuration where a single ziazag 
nanoribbon is connected by two triangular graphene corners, (b) 
Illustration of the Peierls instability of graphene nanoribbons with 
W = 1 and W = 3, where W represents the width. Single and 
double bonds alternatingly appears, (c) The band structure around 
/c = TT of a polyacene with the Peierls instability, (d) The logarithm 
plot of the band gap A, logio A/t, of nanoribbons with various W. 
(e) The logarithm plot of the band gap A, logio A/t, of nanoribbons 
as a function of the applied force S for various W. 

ing an advantage of this property, we propose a nanomechan- 
ical switch sensor. 

The model Hamiltonian is the Su-Schrieffer-Heeger-like 
modeli^, 

^o = ^t,,4aC„a. (19) 

In this model, when the Peierls instability occurs, the trans- 
fer integral takes two values corresponding to the single and 
double bonds, and otherwise it takes one value corresponding 
to the conjugate bond. We are able to calculate the gap ana- 
lytically in the case of polyacene (W = 1), where the band 
structure around k = n is given by 

e (k:,) = ±i (^-t2 + ^4t? + 5^2 + 8^1^2 cos k^^ . (20) 

The band gap is determined by substituting kx = tt sls 



(21) 

We show e (kx) in Fig. 10(c) by taking the values ti = t — 
St I = OMt for single bonds and t2 = t ^ 5t2 = IMt for 
double bonds. We have carried out a numerical estimation 



of the band sturacture for > 2. A tiny gap A opens at 
k = TT. The logarithm plot of the band gap A as a function 
of the width W is shown in Fig. 10(d). The gap decreases 
exponentially as a function of the width, ex 10"^ ''^. The gap 
of the case W = 1 (polyacene) is 49meV, and that of the case 
W = 2is O.SmeV. 

We next compare the energy gain from this gap with the en- 
ergy cost from the elastic energy of a lattice deformation. The 
ground- state energy difference between distorted and undis- 
torted structures is very tiny in polyacene^^. On the other 
hand, the elastic energy cost is proportional to the width W, 
and hence the elastic energy cost becomes larger than the gap 
energy gain for wider nanoribbons. The Peierls instability will 
not occur spontaneously in nanoribbons with W > 2. 

We have also estimated how the band gap A depends on 
the external force. The transfer integrals change proportion- 
ally to the external force. For simplicity we have set ti = 
(1 - 0.06A) t and ^2 = (1 + 0.08A) t. It takes about lOGPa 
for deformation 0.0 It for graphene, while it takes about IGPa 
for deformation 0.0 It for narrow graphene nanoribbons^ ^ We 
show the logarithm plot of the A dependence of the band gaps 
in Fig. 10(e). 

We propose an application of the above system. By strech- 
ing a nanoribbon along the ribbon direction, horizontal bonds 
are stretched and vertical bonds are shrinked. The resultant 
structure is resemble to the deformed structure induced by the 
Peierls transition. We may call it a strain induced Peierls tran- 
sition. When we stretch this structure, the band gap opens and 
the conductance at the zero energy becomes zero by the strain 
induced Peierls-transition. On the other hand, without the ex- 
ternal mechanical force, nanoribbons with > 2 is gapless 
and the system is conductive. We can switch the conductance 
from on to off by stretching the system. The system acts as a 
nanomechanical switch sensor detecting nanoscale displace- 
ment. 



VII. CONCLUSIONS 

A nanodisk can be used as a spin filter just as in a metal- 
ferromagnet-metal junction. A novel feature is that the direc- 
tion of the spin can be controlled by external field or spin cur- 
rent. We have newly proposed nanodisk arrays and nanome- 
chanical switch. These nanodisk-nanoribbon complex struc- 
ture will open a new field of nanoelectronics, spintronics and 
nanoelectromechanics purely based on graphene. 
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